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ABSTRACT: We analyzed the possibility to increase density of field-effect transistors framework a 
single-stage differential-drive CMOS rectifier. Based on this analysis we formulate recommendations to 
obtain the above increasing with increasing of homogeneity of dopant distribution in doped area. We 
also consider a possibility to decrease value of mismatch-induced stress in the considered 
heterostructure. We introduce an analytical approach to analyze mass and heat transport in 
heterostructures during manufacturing of integrated circuits with account mismatch-induced stress. 
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INTRODUCTION 


In the present time several actual problems of 
the solid state electronics (such as increasing of 
performance, reliability and density of elements 
of integrated circuits: diodes, field-effect and 
bipolar transistors) are intensively solving [1- 
6]. To increase the performance of these devices 
it is attracted an interest determination of 
materials with higher values of charge carrier’s 


mobility [7-10]. One way to decrease 
dimensions of elements of integrated circuits is 
manufacturing them in thin film 


heterostructures [3-5,11]. In this case it is 
possible to use inhomogeneity of 
heterostructure and necessary optimization of 
doping of electronic materials [12] and 
development of epitaxial technology to improve 
these materials (including analysis of mismatch 
induced stress) [13-15]. An _ alternative 
approaches to increase dimensions of 
integrated circuits are using of laser and 
microwave types of annealing [16-18]. 


Framework the paper we introduce an 
approach to manufacture field-effect 
transistors. The approach gives a possibility to 
decrease their dimensions with increasing their 
density framework a single-stage differential- 
drive CMOS rectifier. We also consider 
possibility to decrease mismatch-induced stress 
to decrease quantity of defects, generated due 
to the stress. In this paper we consider a 
heterostructure, which consist of a substrate 


and an epitaxial layer (see Fig. 1). We also 
consider a buffer layer between the substrate 
and the epitaxial layer. The epitaxial layer 
includes into itself several sections, which were 
manufactured by using another material. This 
section has been doped by diffusion or ion 
implantation to manufacture the required types 
of conductivity (p or n). These areas became 
sources, drains and gates (see Fig. 1). 


Fig. 1a Structure of the considered rectifier [19] 


After this doping it is required annealing of 
dopant and/or radiation defects. Main aim of 
the present paper is analysis of redistribution of 
dopant and radiation defects to determine 
conditions, which correspond to decreasing of 
elements of the considered rectifier and at the 
same time to increase their density. At the same 
time we consider a possibility to decrease 
mismatch-induced stress. 
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Fig. 1b. Heterostructure with a substrate, epitaxial layers and buffer layer (view from side) 


METHOD OF SOLUTION 

To solve our aim we determine and analyzed 
spatio-temporal distribution of concentration of 
dopant in the considered heterostructure. We 


determine the distribution by solving the 
second Fick's law in the following form [1,20- 
23] 
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Here C(x,y,z,t) is the spatio-temporal 
distribution of concentration of dopant; Q is the 
atomic volume of dopant; V; is the symbol of 


L, 
surficial gradient; (C(x, y, z,t)dz is the 
0 


surficial concentration of dopant on interface 
between layers of heterostructure (in this 
situation we assume, that Z-axis is 
perpendicular to interface between layers of 
heterostructure); tW1(x,y,z,t) and b2(x,y,z,t) are 
the chemical potential due to the presence of 


mismatch-induced stress and _ porosity of 
material; D and Ds are the coefficients of 
volumetric and surficial diffusions. Values of 
dopant diffusions coefficients depends on 
properties of materials of heterostructure, 
speed of heating and cooling of materials during 
annealing and spatio-temporal distribution of 
concentration of dopant. Dependences of 
dopant diffusions coefficients on parameters 
could be approximated by the following 
relations [24-26] 


C"(x, y, Z,t) V (x, y, Zt) V(x, y, Zt) 
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Here D, (x,y,z,T) and Djs (x,y,z,T) are the spatial 
(due to accounting all layers’ of 
heterostruicture) and temperature (due to 
Arrhenius law) dependences of dopant diffusion 
coefficients; T is the temperature of annealing; P 
(xy,z,T) is the limit of solubility of dopant; 
parameter y depends on properties of materials 
and could be integer in the following interval y 
€[1,3] [24]; V (xy,z,t) is the spatio-temporal 
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PY xyz.) 
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distribution of concentration of radiation 
vacancies; V* is the equilibrium distribution of 
vacancies. Concentrational dependence of 
dopant diffusion coefficient has been described 
in details in [24]. Spatio-temporal distributions 
of concentration of point radiation defects have 
been determined by solving the following 
system of equations [20-23,25,26] 
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Here I! (x,y,z,t) is the spatio-temporal 
distribution of concentration of radiation 
interstitials; /* is the equilibrium distribution of 
interstitials; D)(x,y,z,T), Dv(x,y,z,T), Dis(x,y, Z,T), 
Dys(Xy,Z,T) are the coefficients of volumetric 
and surficial diffusions of interstitials and 
vacancies, respectively; terms V2(x,y,z,t) and 
2(xy,z,t) correspond to generation § of 
divacancies and di interstitials, respectively 
(see, for example, [26] and appropriate 
references in this book); kry(x,y,z,T), Kri(x,y,zZ,T) 
and kyy(xy,z,T) are the parameters’ of 
recombination of point radiation defects and 
generation of their complexes; k is the 
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Boltzmann constant; @ = a3, a is the interatomic 


distance; ¢ is the specific surface energy. To 
account porosity of buffer layers we assume, 
that porous are approximately cylindrical with 


2 2 
average values F=./X, + y, and Z: before 


annealing [23]. With time small pores are 
decomposing on vacancies. These vacancies are 
absorbing by larger pores [27]. With time large 
pores became larger due to absorbing the 
vacancies and became more spherical [27]. 
Distribution of concentration of vacancies in 
heterostructure, existing due to porosity, could 
be determined by summing on all pores, i.e. 


lm on 
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Here a, B and x are the average distances 
between centers of pores in directions x, y and 
z; Il, m and n are the quantity of pores 
inappropriate directions. Spatio-temporal 


distributions of divacancies ®y (x,,z,t) and di 
interstitials ®; (x,y,z, t) could be determined by 
solving the following system of equations 
[25,26] 
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Here Doi(x,y,z,T), Dov(Xy,Z,T), Deis (X,y,Z,T) and Devs(x,y,z,T) are the coefficients of volumetric and 
surficial diffusions of complexes of radiation defects; k;(x,y,z,T) and ky(x,y,z,T) are the parameters of 
decay of complexes of radiation defects. 


Chemical potential 111 in Eq.(1) could be determine by the following relation [20] 
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tensor; uj uj; are the components u;(x,y,z,t), Uy(xy,z,t) and u,(x,,z,t) of the displacement vector 
Uu (x, y, Zt ); X;, Xj are the coordinate x, y, z. The Eq. (3) could be transform to the following form 
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is the mismatch parameter; as, dz, are lattice (for our case) with room _ temperature. 
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B is the coefficient of thermal expansion; T, is [21] 
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Conditions for the system of Eq. (8) could be written in the form 
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We determine spatio-temporal distributions of concentrations of dopant and radiati-on defects by 
solving the Eqs.(1), (3) and (5) framework standard method of averaging of function corrections [28]. 


Previously we transform the Eqs.(1), (3) and (5) to the following form with account initial 
distributions of the considered concentrations 
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order approximations of the required concentrations in the following form 
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Integration of the left and right sides of the Eqs. (1b), (3b) and (5b) on time gives us possibility to 
obtain relations for above approximation in the final form 
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We determine average values of the first-order approximations of concentrations of dopant and 
radiation defects by the following standard relation [28] 
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Integration of the left and the right sides of Eqs. (1d), (3d) and (5d) gives us possibility to obtain 
relations for the required concentrations in the final form 
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Average values of the second-order approximations of required approximations by using the following 
standard relation [28] 
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Approximations of the second and higher orders of components of displacement vector could be 
determined by using standard replacement of the required components on the following sums 
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Integration of the left and right sides of the above relations on time t leads to the following result 


u SE(:) \2 a [Ju (x y,Z, t)dt t+ 15] R()- 


uj: y.24)= Syke 6[1+0 (z)|J A 


E(z) a2 18 E(z) Bias 
9 9 9 d d& : : : d de 
Sica Gy aay ivtovedaeadors TA sunlinac)acao. 


o flu(x»20)dtao | 4 © thy, (ny,2n)ardo 1K (2) 
tage hie nz . aeO) sOe.” iia cai al 


E(z) B (z) O18 Q? ~6 
Tal K(2)E aa) da sss emer at dd 7 [Jal %t dt dd x 
1 E(z) Ae ose 
Ke 6[l+o ee [KO-5 oe pao J Ju, (x y.zt)dt dd x 


Pankratov 


www.mahendrapublications.com 


International Journal of Advanced Science and Engineering 


Int. J. Adv. Sci. Eng. Vol.8 No.4. 2391-2417 (2022) 2411 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


2 ee. 2 
x ag Sl Y,Z,0 t)de d+ [lu (x Y, ZU s)ae do 


x 


Zo) 


= +2) cA ie X,Y, 2,0 )dt dB +uUy, + B (z) 
KG) false! J mI ys )d dd + Uo, K(z) (z) 


oot 
x" Fir (x,y,2z0) de dd 
0x00 


2 19 td 
uy,(x, y, Zt) = ls, o[fu.(.y,20)dtdd+— fun yet )dt do | x 
Ox’ 00 OxO yoo 


2 
1 K(z Oz 18 1 5E(z) 
x + JJu,,(% y.2.t dt 10+ pee + Kos 
2 ce 


| fie, (x, 9,20) atddo+—— a E(z) 2 [u(x y.2t) dtd0+ 
) 


on 00 2p (z) Az 1+ o (z) OZ00 


6) B(z to E(z) 

+s LY, 2,0 dt || Ee) JIT (x, y,z,t dt #o-| ee 
1 td E(z) a? ae 

-K( sacagl lino y> z,t )dt dvd eth = J fun (2, y, 2,0 )dt aeve+ 


il L__x(B@fir _KQ), 
+52 fhlssner)ardo] ae a ,y.z,t)dtdd 0) 


ral ov 1 ral Pan 5E(z) 
AVA »y,z,t )dt dd ——_._ —— ¥,z,0) dtav< —— 
Bxdy iit y z,t) T p(z)day’ jfu.(x y z,t) T | ape 
O | Elz) | oO Ae 
. »y,Z,t) dt dB +— ,y,z,t)dtdd 
° a} Aol eline Yast) dt rte y,z,t) dt } 
1 1 E(z) a? i 
ele yZt) dt dd +uy, 
"2p (z) nah (:) eat aas lint y ae) G +Up, 


Elz CG? =e G* «0 
u(x, y, Z1)= eal = JJu.(x, 2%) dtd8 + By feels y,Z,T) dtd&+ 


Pankratov 


International Journal of Advanced Science and Engineering www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.8 No.4. 2391-2417 (2022) 2412 


E-ISSN: 2349 5359; P-ISSN: 2454-9967 


avd 


{ Ju, (x, y> en) 
00 


dtd&d cg : 


p(z) plz) 


x 


ov ov 
g [ fu, (x, y, 2,7) dt tie a y,z,0) dtd&+ 
00 Oyoo 


E(z) 6 0 


a. ie ' 
Fo MAX ZT dtd+ 
ae nt Y,Z ) Bayz 
O 
x —< K as 
cS (2 
ot 
+—fJu,,(xy,z,t) dtd ee | o 
OZ 00 


oot 
- Fh. y, z,T) dtdv- 


X00 


Cc 


B(z) é 


_K ( zy 
p (z) OZ 00 
Framework this paper we _ determine 
concentration of dopant, concentrations of 
radiation defects and components _ of 


displacement vector by using the second-order 
approximation framework method of averaging 
of function corrections. This approximation is 
usually enough good approximation to make 
qualitative analysis and to obtain some 
quantitative results. All obtained results have 
been checked by comparison with results of 
numerical simulations. 


DISCUSSION 

In this section we analyzed dynamics of 
redistributions of dopant and radiation defects 
during annealing and under influence of 
mismatch-induced stress and modification of 
porosity. Typical distributions of concentrations 
of dopant in heterostructures are presented on 
Figs. 2 and 3 for diffusion and ion types of 
doping, respectively. These distributions have 
been calculated for the case, when value of 
dopant diffusion coefficient in doped area is 
larger, than in nearest areas. The figures show, 


0.54 
Epitaxial laver 


6p (z) Oz [1+6 (z) 


[fay y,z,t)dtdb— 
OZ00 


ra) ot ra) ov 
— | Ju,,(x, y,z,t) dtd8—-— Jfu,,(x, y,z,t) dt dd - 
Oy 00 0Z00 


HIT (xy,et) dt d9 + Uy, 


that inhomogeneity of heterostructure gives us 
possibility to increase compactness’ of 
concentrations of dopants and at the same time 
to increase homogeneity of dopant distribution 
in doped part of epitaxial layer. However 
framework this approach of manufacturing of 
bipolar transistor it is necessary to optimize 
annealing of dopant and/or radiation defects. 
Reason of this optimization is following. If 
annealing time is small, the dopant did not 
achieve any interfaces between materials of 
heterostructure. In this situation one cannot 
find any modifications of distribution of 
concentration of dopant. If annealing time is 
large, distribution of concentration of dopant is 
too homogenous. We optimize annealing time 
framework recently introduces approach [29- 
39]. Framework this criterion we approximate 
real distribution of concentration of dopant by 
step-wise function (see Figs. 4 and 5). Farther 
we determine optimal values of annealing time 
by minimization of the following mean-squared 


error 


Substrate 
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Fig.2. Distributions of concentration of infused dopant in heterostructure from Fig. 1 in direction, which is 
perpendicular to interface between epitaxial layer substrate. Increasing of number of curve corresponds to 
increasing of difference between values of dopant diffusion coefficient in layers of heterostructure under 
condition, when value of dopant diffusion coefficient in epitaxial layer is larger, than value of dopant diffusion 
coefficient in substrate 
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Fig.3. Distributions of concentration of implanted dopant in heterostructure from Fig. 1 in direction, which is 
perpendicular to interface between epitaxial layer substrate. Curves 1 and 3 corresponds to annealing time © = 
0.004.8(Lx2+Ly2+Lz2)/Do. Curves 2 and 4 corresponds to annealing time © = 0.0057 (Lx2+Ly2+Lz2)/Do. Curves 1 and 
2 corresponds to homogenous sample. Curves 3 and 4 corresponds to heterostructure under condition, when 


value of dopant diffusion coefficient in epitaxial layer is larger, than value of dopant diffusion coefficient in 
substrate 
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Fig. 4. Spatial distributions of dopant in heterostructure after dopant infusion. 


Curve 1 is idealized distribution of dopant. Curves 2-4 are real distributions of dopant for different 
values of annealing time. Increasing of number of curve corresponds to increasing of annealing time 


1 Ly LyL. 


"ELL. rae (x, y,z,0)-y (x, y,z)|dzd ydx, (15) 


where y (x,y,z) is the approximation function. dopant achieves appropriate interfaces between 


Dependences of optimal values of annealing 
time on parameters are presented on Figs. 6 and 
7 for diffusion and ion types of doping, 
respectively. It should be noted, that it is 
necessary to anneal radiation defects after ion 
implantation. One could find spreading of 
concentration of distribution of dopant during 
this annealing. In the ideal case distribution of 


materials of heterostructure during annealing of 
radiation defects. If dopant did not achieve any 
interfaces during annealing of radiation defects, 
it is practicably to additionally anneal the 
dopant. In this situation optimal value of 
additional annealing time of implanted dopant 
is smaller, than annealing time of infused 
dopant 
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C(x,0) 


Fig. 5. Spatial distributions of dopant in heterostructure after ion implantation. Curve 1 is idealized 
distribution of dopant. Curves 2-4 are real distributions of dopant for different values of annealing time. 
Increasing of number of curve corresponds to increasing of annealing time 


O55)» 
0.4 j 
i 
«. 02 
al e 
2 
st 
2 02 
1 
0.1 


0.0 0.1 0.2 0.3 0.4 0.5 
a/L, §, &, Y 


Fig.6. Dependences of dimensionless optimal annealing time for doping by diffusion, which have been obtained 
by minimization of mean-squared error, on several parameters. Curve 1 is the dependence of dimensionless 
optimal annealing time on the relation a/L and € = y = 0 for equal to each other values of dopant diffusion 
coefficient in all parts of heterostructure. Curve 2 is the dependence of dimensionless optimal annealing time 
on value of parameter ¢ for a/L=1/2 and & = y = 0. Curve 3 is the dependence of dimensionless optimal 
annealing time on value of parameter € for a/L=1/2 and € = y = 0. Curve 4 is the dependence of dimensionless 
optimal annealing time on value of parameter y for a/L=1/2 and ec =&=0 
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Fig.7. Dependences of dimensionless optimal annealing time for doping by ion implantation, which have been 
obtained by minimization of mean-squared error, on several parameters. Curve 1 is the dependence of 
dimensionless optimal annealing time on the relation a/L and € = y = 0 for equal to each other values of dopant 
diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of dimensionless optimal 
annealing time on value of parameter ¢ for a/L=1/2 and € = y = 0. Curve 3 is the dependence of dimensionless 
optimal annealing time on value of parameter & for a/L=1/2 and € = y = 0. Curve 4 is the dependence of 
dimensionless optimal annealing time on value of parameter y for a/L=1/2 and ¢ = =0 
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Fig. 8. Normalized dependences of component u, of displacement vector on coordinate z for nonporous 


(curve 1) and porous (curve 2) epitaxial layers 


Further we analyzed influence of relaxation of 
mechanical stress on distribution of dopant in 
doped areas of heterostructure. Under following 
condition €9< O one can find compression of 
distribution of concentration of dopant near 
interface between materials of heterostructure. 
Contrary (at €0>0) one can find spreading of 
distribution of concentration of dopant in this 
area. This changing of distribution of 
concentration of dopant could be at least 
partially compensated by using laser annealing 
[37]. This type of annealing gives us possibility 
to accelerate diffusion of dopant and another 
process in annealed area due to inhomogenous 
distribution of temperature and Arrhenius law. 

0.4 


V(z) 0.2 4 


0.0 


Accounting relaxation of mismatch-induced 
stress in heterostructure could leads _ to 
changing of optimal values of annealing time. At 
the same time modification of porosity gives us 
possibility to decrease value of mechanical 
stress. On the one hand mismatch-induced 
stress could be used to increase density of 
elements of integrated circuits. On the other 
hand could leads to generation dislocations of 
the discrepancy. Figs. 8 and 9 show 
distributions of concentration of vacancies in 
porous materials and component of 
displacement vector, which is perpendicular to 
interface between layers of heterostructure. 


0.0 


1 
a 
z 


Fig. 9. Normalized dependences of vacancy concentrations on coordinate z in unstressed (curve 1) 
and stressed (curve 2) epitaxial layers 


CONCLUSION 

In this paper we model redistribution of 
infused and implanted dopants with account 
relaxation mismatch-induced stress during 
manufacturing field-effect heterotransistors 
framework a_ single-stage differential-drive 
CMOS rectifier. We formulate recommendations 
for optimization of annealing to decrease 
dimensions of transistors and to increase their 
density. We formulate recommendations to 
decrease mismatch-induced stress. Analytical 
approach to model diffusion and ion types of 


doping with account concurrent changing of 
parameters in space and time has been 
introduced. At the same time the approach gives 
us possibility to take into account nonlinearity 
of considered processes. 
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